BAR4R-PUB-10/017 
SLAC-PUB-14270 



Search for the Decay B° — > 77 

P. del Amo Sanchez, J. P. Lees, V. Poireau, E. Prencipe, and V. Tisserand 
Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 

J. Garra Tico and E. Grauges 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

M. Martinelli ab , A. Palano ab , and M. Pappagallo a& 
INFN Sezione di Bari a ; Dipartimento di Fisica, Universita di Bari b , 1-10126 Bari, Italy 

G. Eigen, B. Stugu, and L. Sun 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

M. Battaglia, D. N. Brown, B. Hoobcrman, L. T. Kerth, 
Yu. G. Kolomensky, G. Lynch, I. L. Osipenkov, and T. Tanabe 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

CM. Hawkes and A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

H. Koch and T. Schroeder 
Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 

D. J. Asgeirsson, C. Hearty, T. S. Mattison, and J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 

A. Khan and A. Randle-Conde 

Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

V. E. Blinov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, A. P. Onuchin, 
S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, K. Yu. Todyshev, and A. N. Yushkov 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioli, S. Curry, D. Kirkby, A. J. Lankford, M. Mandelkern, E. C. Martin, and D. P. Stoker 
University of California at Irvine, Irvine, California 92697, USA 

H. Atmacan, J. W. Gary, F. Liu, O. Long, and G. M. Vitug 

University of California at Riverside, Riverside, California 92521, USA 

C. Campagnari, T. M. Hong, D. Kovalskyi, J. D. Richman, and C. West 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

A. M. Eisner, C. A. Heusch, J. Kroseberg, W. S. Lockman, 
A. J. Martinez, T. Schalk, B. A. Schumm, A. Seiden, and L. O. Winstrom 
University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

C. H. Cheng, D. A. Doll, B. Echenard, D. G. Hitlin, P. Ongmongkolkul, F. C. Porter, and A. Y. Rakitin 

California Institute of Technology, Pasadena, California 91125, USA 

R. Andreassen, M. S. Dubrovin, G. Mancinelli, B. T. Meadows, and M. D. Sokoloff 

University of Cincinnati, Cincinnati, Ohio 45221, USA 



2 



P. C. Bloom, W. T. Ford, A. Gaz, M. Nagel, U. Nauenberg, J. G. Smith, and S. R. Wagner 
University of Colorado, Boulder, Colorado 80309, USA 

R. Ayad3 and W. H. Toki 

Colorado State University, Fort Collins, Colorado 80523, USA 

H. Jasper, T. M. Karbach, J. Merkel, A. Pctzold, B. Spaan, and K. Wacker 
Techmsche Universitat Dortmund, Fakultat Physik, D-44221 Dortmund, Germany 

M. J. Kobel, K. R. Schubert, and R. Schwierz 
Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard and M. Verderi 
Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Poly technique, F-91128 Palaiseau, France 

P. J. Clark, S. Playfer, and J. E. Watson 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andrcotti ab , D. Bettoni a , C. Bozzi a , R. Calabrese af \ A. Cecchi ab , G. Cibinetto ab , E. Fioravanti ab , 
P. Franchini ab , E. Luppi ab , M. Munerato ab , M. Negrini" 6 , A. Petrella ab , and L. Piemontese a 
INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universita di Ferrara b , 1-44100 Ferrara, Italy 

R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, M. Nicolaci, 
S. Pacetti, P. Patteri, I. M. PeruzziQ M. Piccolo, M. Rama, and A. Zallo 
INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

R. Contri ab , E. Guido ab , M. Lo Vetcre ab , M. R. Monge ab , S. Passaggio Q , C. Patrignani Qb , E. Robutti Q , and S. Tosi a ' 
INFN Sezione di Genova" '; Dipartimento di Fisica, Universita di Genova b , 1-16146 Genova, Italy 

B. Bhuyan and V. Prasad 
Indian Institute of Technology Guwahati, Guwahati, Assam, 781 039, India 

C. L. Lee and M. Morii 

Harvard University, Cambridge, Massachusetts 02138, USA 

A. Adametz, J. Marks, and U. Uwer 
Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

F. U. Bernlochner, M. Ebert, H. M. Lacker, T. Lueck, and A. Volk 
Humboldt-Universitat zu Berlin, Institut fur Physik, Newtonstr. 15, D-12489 Berlin, Germany 

P. D. Dauncey and M. Tibbetts 
Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Behera and U. Mallik 

University of Iowa, Iowa City, Iowa 52242, USA 

C. Chen, J. Cochran, H. B. Crawley, L. Dong, W. T. Meyer, S. Prell, E. I. Rosenberg, and A. E. Rubin 

Iowa State University, Ames, Iowa 5001 1-3160, USA 

A. V. Gritsan and Z. J. Guo 
Johns Hopkins University, Baltimore, Maryland 21218, USA 

N. Arnaud, M. Davier, D. Derkach, J. Firmino da Costa, G. Grosdidier, F. Le Diberder, A. M. Lutz, B. Malaescu, 
A. Perez, P. Roudeau, M. H. Schune, J. Serrano, V. Sordini^ A. Stocchi, L. Wang, and G. Wormser 
Laboratoire de V Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 



3 



D. J. Lange and D. M. Wright 

Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

I. Bingham, C. A. Chavez, J. P. Coleman, J. R. Fry, E. Gabathuler, 
R. Garnet, D. E. Hutchcroft, D. J. Payne, and C. Touramanis 
University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bcvan, F. Di Lodovico, R. Sacco, and M. Sigamani 
Queen Mary, University of London, London, El 4NS, United Kingdom 

G. Cowan, S. Paramesvaran, and A. C. Wren 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. N. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

A. G. Denig, M. Fritsch, W. Gradl, and A. Hafner 

Johannes Gutenberg-Universitdt Mainz, Institut fur Kernphysik, D-55099 Mainz, Germany 

K. E. Alwyn, D. Bailey, R. J. Barlow, G. Jackson, and G. D. Laffcrty 
University of Manchester, Manchester M13 9PL, United Kingdom 

J. Anderson, R. Cenci, A. Jawahery, D. A. Roberts, G. Simi, and J. M. Tuggle 

University of Maryland, College Park, Maryland 20742, USA 

C. Dallapiccola and E. Salvati 

University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, D. Dujmic, G. Sciolla, and M. Zhao 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

D. Lindemann, P. M. Patel, S. H. Robertson, and M. Schram 
McGill University, Montreal, Quebec, Canada H3A 2T8 

P. Biassoni afc , A. Lazzaro afc , V. Lombardo a , F. Palombo afc , and S. Stracka ab 
INFN Sezione di Milano a ; Dipartimento di Fisica, Universitd di Milano b , 1-20133 Milano, Italy 

L. Cremaldi, R. GodangH R. Kroeger, P. Sonnek, and D. J. Summers 
University of Mississippi, University, Mississippi 38677, USA 

X. Nguyen, M. Simard, and P. Taras 
Universitd de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 

G. De Nardo afc , D. Monorchio afc , G. Onorato ah , and C. Sciacca afc 
INFN Sezione di Napoli a ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico IF, 1-80126 Napoli, Italy 

G. Raven and H. L. Snoek 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

C. P. Jessop, K. J. Knoepfel, J. M. LoSecco, and W. F. Wang 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

L. A. Corwin, K. Honscheid, R. Kass, and J. P. Morris 
Ohio State University, Columbus, Ohio 43210, USA 

N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, R. Rahmat, N. B. Sinev, D. Strom, J. Strubc, and E. Torrence 

University of Oregon, Eugene, Oregon 97403, USA 



4 



G. Castelli ab , E. Feltresi ab , N. Gagliardi ab , M. Margoni ab , M. Morandin , 
M. Posocco Q , M. Rotondo a , F. Simonetto af \ and R. Stroili Qb 
INFN Sezione di Padova a ; Dipartvmento di Fisica, Universitd di Padova b , 1-35131 Padova, Italy 

E. Ben-Haim, G. R. Bonneaud, H. Briand, G. Calderini, J. Chauveau, 
O. Hamon, Ph. Lcrustc, G. Marchiori, J. Ocariz, J. Prendki, and S. Sitt 

Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Mane Curie-Paris6, 
Universite Denis Diderot-Paris'} ', F-75252 Paris, France 

M. Biasini ab , E. Manoni ab , and A. Rossi ab 
INFN Sezione di Perugia a ; Dipartimento di Fisica, Universitd, di Perugia b , 1-06100 Perugia, Italy 

C. Angelini ab , G. Batignani ab , S. Bettarini ab , M. Carpinelli ab G. Casarosa ab , A. Cervelli ab , 
F. Forti ab , M. A. Giorgi ab , A. Lusiani ac , N. Neri ab , E. Paoloni ab , G. Rizzo ab , and J. J. Walsh a 
INFN Sezione di Pisa a ; Dipartimento di Fisica, Universitd di Pisa b ; Scuola Norrnale Superiore di Pisa c , 1-56127 Pisa, Italy 

D. Lopes Pegna, C. Lu, J. Olsen, A. J. S. Smith, and A. V. Telnov 

Princeton University, Princeton, New Jersey 08544, USA 

F. Anulli a , E. Baracchini ab , G. Cavoto a , R. Faccini ab , F. Ferrarotto a , F. Ferroni ab , 

M. Gaspero Qb , L. Li Gioi a , M. A. Mazzoni a , G. Piredda a , and F. Renga Qb 
INFN Sezione di Roma a ; Dipartimento di Fisica, 
Universitd di Roma La Sapienza b , 1-00185 Roma, Italy 

T. Hartmann, T. Leddig, H. Schroder, and R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, B. Franek, E. O. Olaiya, and F. F. Wilson 

Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, G. Hamel de Monchenault, G. Vasseur, Ch. Yeche, and M. Zito 
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 

M. T. Allen, D. Aston, D. J. Bard, R. Bartoldus, J. F. Benitez, C. Cartaro, M. R. Convery, J. Dorfan, 
G. P. Dubois-Felsmann, W. Dunwoodie, R. C. Field, M. Franco Sevilla, B. G. Fulsom, A. M. Gabareen, 

M. T. Graham, P. Grenier, C. Hast, W. R. Innes, M. H. Kelsey, H. Kim, P. Kim, M. L. Kocian, 
D. W. G. S. Leith, S. Li, B. Lindquist, S. Luitz, V. Luth, H. L. Lynch, D. B. MacFarlane, H. Marsiske, 
D. R. Muller, H. Neal, S. Nelson, C. P. O'Grady, I. Ofte, M. Perl, T. Pulliam, B. N. Ratcliff, 
A. Roodman, A. A. Salnikov, V. Santoro, R. H. Schindler, J. Schwiening, A. Snyder, D. Su, 
M. K. Sullivan, S. Sun, K. Suzuki, J. M. Thompson, J. Va'vra, A. P. Wagner, M. Weaver, C. A. West, 
W. J. Wisniewski, M. Wittgen, D. H. Wright, H. W. Wulsin, A. K. Yarritu, C. C. Young, and V. Ziegler 
SLAC National Accelerator Laboratory, Stanford, California 94309 USA 

X. R. Chen, W. Park, M. V. Purohit, R. M. White, and J. R. Wilson 

University of South Carolina, Columbia, South Carolina 29208, USA 

S. J. Sekula 

Southern Methodist University, Dallas, Texas 75275, USA 

M. Bellis, P. R. Burchat, A. J. Edwards, and T. S. Miyashita 
Stanford University, Stanford, California 94305-4060, USA 

S. Ahmed, M. S. Alam, J. A. Ernst, B. Pan, M. A. Saeed, and S. B. Zain 

State University of New York, Albany, New York 12222, USA 



N. Guttman and A. Soffer 



5 



Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 

P. Lund and S. M. Spanier 
University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. M. Ruland, C. J. Schilling, R. F. Schwitters, and B. C. Wray 
University of Texas at Austin, Austin, Texas 78712, USA 

J. M. Izen and X. C. Lou 

University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi ah , D. Gamba ab , and M. Pelliccioni afc 
INFN Sezione di Torino 0, ; Dipartimento di Fisica Sperimentale, Universita di Torino" , 1-10125 Torino, Italy 

M. Bomben afc , L. Lanceri Q& , and L. Vitale Q& 
INFN Sezione di Trieste"" ; Dipartimento di Fisica, Universita di Trieste b , 1-34-127 Trieste, Italy 

N. Lopez-March, F. Martinez- Vidal, D. A. Milanes, and A. Oyanguren 
IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spain 

J. Albert, Sw. Banerjee, H. H. F. Choi, K. Hamano, G. J. King, 
R. Kowalewski, M. J. Lewczuk, I. M. Nugent, J. M. Roney, and R. J. Sobie 
University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

T. J. Gershon, P. F. Harrison, T. E. Latham, and E. M. T. Puccio 

Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, S. Dasu, K. T. Flood, Y. Pan, R. Prepost, C. O. Vuosalo, and S. L. Wu 

University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: October 13, 2010) 

We report the result of a search for the rare decay B° — ► 77 in 426 fb _1 of data, corresponding 
to 226 million B°B° pairs, collected on the T(45) resonance at the PEP-II asymmetric-energy e e 
collider using the BABAR detector. We use a maximum likelihood fit to extract the signal yield and 
observe 211JJ signal events with a statistical signficance of 1.9 a. This corresponds to a branching 
fraction B(B° ->■ 77) = (1.7 ± l.l(stat.) ± 0.2(syst.)) x 10" 7 . Based on this result, we set a 90% 
confidence level upper limit of B{B° ->■ 77) < 3.3 x 10~ 7 . 

PACS numbers: 13.20.He 



I. INTRODUCTION 

In the standard model (SM), the decay B° — » 77 oc- 
curs through a flavor-changing neutral current (FCNC) 
transition involving electroweak loop diagrams, as illus- 
trated in Figure [TJ The decays B° —> 77 and B s —> 77 
are closely related, with the b — ¥ ^77 transition be- 
ing suppressed with respect to b — > S77 by Cabbibo- 
Kobayashi-Maskawa (CKM) factors {\V td \ 2 /\V ts \ 2 ~ 
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0.04). Hadron dynamics introduces uncertainties into 
the prediction of branching fractions for these decays and 
may modify the ratio away from the CKM-implied value. 
While B° — » 77 is expected to have a smaller branching 
fraction than B s — > 77, a search for the latter faces the 
experimental challenge of obtaining a large sample of B s 
mesons, whereas large samples of B° mesons are read- 
ily available from B Factory experiments running on the 
T(iS) resonance. 

A leading order calculation for the branching fraction 
of B° 77 yields an estimate of 3.1+^6 x 10~ 8 1]. This 
mode is sensitive to new physics that could lead to an 
enhancement of the branching fraction due to possible 
contributions of non-SM heavy particles occurring in the 
loop of the leading-order Feynman diagrams. Such en- 
hancements to the branching fraction for B° — » 77 are 
less constrained than those for B s — > 77 due to the fact 
that the b — > sj transition, responsible for B s — > 77, is 
known much more accurately than b —¥ dr/. For example, 
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FIG. 1: Examples of lowest order SM Feynman diagrams for 
B° — > 77. The symbol q represents a u, c, or t quark. In 
some new physics scenarios, the VK-boson may be replaced by 
a charged Higgs particle. 



some new physics scenarios involving an extended Higgs 
sector may considerably enhance the branching fractions 
with respect to the SM expectation 0]. Supersymmetry 
with broken i?-parity Q also provides scenarios where or- 
der of magnitude enhancements are possible. In addition, 
since the two-photon final state can be either CP-even or 
CP-odd, studies of CP-violating effects may ultimately 
be possible. 

The best previous upper limit on the branching frac- 
tion at 90% confidence level (CL) is B(B° ->■ 77) < 
6.2 x 10~ 7 set by the Belle experiment [H using a dataset 
recorded at the T(45) resonance with an integrated lu- 
minosity of 104 fb _1 . For the related process B s — > 77, 
Belle has set a an upper limit on the branching fraction 
of 8.7 x 10~ 6 (90% CL) based on 23.6 ftr 1 of data taken 
on the T(5S) resonance 

We report herein a new search for the decay B° — > 77 
which uses a data sample with integrated luminosity of 
426 ftV 1 taken at the T(AS) resonance. This corresponds 
to the entire BABAR T(4S) dataset and contains 226 mil- 
lion B°B° pairs. The analysis does not distinguish be- 
tween B° and B°, and throughout this article, charge 
conjugation is implied for all reactions. 

The analysis proceeds through several steps. The full 
dataset is first reduced to a manageable size by selecting 
events based on loose kinematic criteria consistent with 
the B° — > 77 hypothesis. Studies are then performed 
to determine an optimal set of event selection criteria to 
maximize the efficiency for detecting B° — > 77 events 
while effectively rejecting background events. This anal- 
ysis is performed "blind" in the sense that the event selec- 
tion criteria are determined without considering the on- 
resonance data within a specific signal region, as denned 
below. We use an unbinned extended maximum like- 
lihood fit to extract the signal yield from the remaining 
events. Finally, since no statistically significant B° — > 77 
signal is observed, an upper limit on the branching frac- 
tion is calculated based on the likelihood function deter- 
mined from the fit. 

Section II of this article describes the BABAR detec- 
tor and the T(AS) dataset. Section III outlines the op- 
timization of the event selection criteria and discusses 
backgrounds due to exclusive B decays. Section IV de- 
scribes the fit methodology and Section V discusses the 
sources of systematic uncertainties. Finally, Section VI 
reports the resulting branching fraction and the upper 



limit. 



II. THE BABAR DETECTOR AND DATASET 

The BABAR detector recorded data from the PEP-II B 
Factory located at the SLAC National Accelerator Labo- 
ratory. In PEP-II, head-on collisions of 9.0 GeV elec- 
trons and 3.1 GeV positrons provide a center-of-mass 
(CM) energy of 10.58 GeV that lies at the peak of the 
T(4S) resonance. The T(4S) meson decays almost ex- 
clusively to BB pairs. The subsequent B meson decays 
were observed in the BABAR detector, which has been de- 
scribed in detail elsewhere [6j . Briefly, a superconducting 
solenoid produces a 1.5 T magnetic field approximately 
parallel to the colliding electron and positron beams. 
Most detector subsystems are inside the volume of the 
solenoid. Charged particle tracking is accomplished by a 
combination of a five layer double-sided silicon strip ver- 
tex detector and a 40 layer drift chamber. The track 
reconstruction algorithm accepts tracks with a trans- 
verse momentum greater than 50 MeV/c. Identification 
of charged particles is accomplished using a ring-imaging 
Cherenkov detector augmented with energy loss measure- 
ments in the tracking detectors. Photons are detected in 
the electromagnetic calorimeter (EMC), which consists 
of 6580 CsI(Tl) crystals, oriented in a quasi-projective 
geometry with respect to the interaction point. Out- 
side of the solenoid, the steel flux return is instrumented 
with a combination of resistive-plate chambers and lim- 
ited streamer tubes to provide detection of muons and 
long-lived neutral hadrons. 

The EMC is the most important detector subsystem 
for the B° — > 77 search. It provides polar angle cov- 
erage in the laboratory frame from 15.8° to 141.8°, and 
full azimuthal coverage, corresponding to a solid-angle 
coverage of 90% in the T(45 l ) CM frame. When a pho- 
ton interacts with the EMC it creates an electromagnetic 
shower depositing its energy into many contiguous crys- 
tals, typically 10 to 30, hereafter called a "cluster" . If no 
track in the event points to the cluster, it is designated 
as a photon candidate. Individual crystals are read out 
by a pair of silicon PIN photodiodes attached to the rear 
face of the crystal. Amplified signals are sampled and 
digitized with a period of 270 ns, providing a continu- 
ous data stream to which gain and pedestal corrections 
are applied. When a first level trigger is recorded, data 
samples in a time window of ±1 /is are selected, produc- 
ing a waveform which is analyzed by a feature extraction 
algorithm running in real time in readout modules. For 
events passing a higher level trigger, any EMC signal with 
energy above a 0.5 MeV threshold has its deposited en- 
ergy and timing information recorded for offline analysis. 

The energy resolution of the EMC is parameterized as 
the sum of two terms added in quadrature, given by Q: 



a E _ 2.30% 
~~E ~ $JE ( GeV) 



8 1.35%, 
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while the angular resolution in the polar angle 9 and the 
azimuthal angle <f> is given by: 



4.16 



■. mrad. 



In addition to the on-resonance data, a data sample 
of 44 ftT 1 taken about 40 MeV below the T(45) peak 
is recorded and used to validate the Monte Carlo (MC) 
simulation of continuum processes, e + e~ — > qq, where 
q = u,d,s, or c) and e + e~ — > t + t~. MC simulated 
events are produced using the EVTGEN Q package to 
model the physics of B meson decays and JETSETj9| to 
model quark fragmentation. The GEANT4 toolkit [l0( is 
used to simulate the interaction of these events with the 
detector model. These tools are designed to take into ac- 
count the varying detector and beam conditions encoun- 
tered during data-taking. The MC events were analyzed 
with the same reconstruction algorithms, event selection 
and fitting procedures as data. MC samples of BB events 
correspond to about four times the integrated luminos- 
ity of the data. Those for e + e~ — > cc events correspond 
to twice the data luminosity, while those for e + e~ — > qq 
(where q = u, d, ors) and e + e~ — > t + t~ correspond to 
approximately the same luminosity as data. In addition, 
special MC data sets are created in which large samples 
of rare B meson decays are generated for the purpose of 
investigating the signal B° — > 77 decay as well as possible 
backgrounds due to other B decays. 



III. EVENT SELECTION AND BACKGROUNDS 
A. Event Selection 

The full T(45) on-resonance dataset is first reduced by 
selecting events that contain at least two photons with 
energies of 1.15 < E* < 3.50 GeV, where the asterisk 
indicates a quantity in the Y(AS) CM frame. We con- 
sider all combinations of two photons whose energies lie 
in this range and add their four-momentum to create B 
meson candidates. Hereafter, these photons are referred 
to as B candidate photons. The distribution of correctly 
reconstructed B candidates will peak in two nearly un- 
corrected variables, toes and AE. The beam energy 
substituted mass is defined as toes = \J £^ am — P§ /c 2 
and the energy difference is AE = E^ — E^ cam , where 
^beam * s the beam energy, and and £3 are the three- 
momentum and energy of the B candidate, respectively. 
For B — > 77 events, the Toes distribution will peak at 
the B meson mass, 5.279 GeV/c 2 [llj]. The MC predicts 
a full- width at half maximum (FWHM) of 6.5MeV/c 2 . 
The AE distribution is asymmetric and will peak near 
zero with a tail to the negative AE side due to photon 
energy loss outside the active volume of the EMC. The 
FWHM for AE predicted by the MC is about 150 MeV. 
We select an event for further analysis if it contains 
exactly one B candidate with toes > 5.1 GeV/c 2 and 



-0.50 < AE < 0.50 GeV. We find that in B° 77 
MC only 0.06% of events have more than one B can- 
didate and are removed by this selection. Events from 
continuum processes, e + e~ — > qq and e + e~ — > t + t~, 
are suppressed by requiring the ratio of the second to 
zeroth Fox- Wolfram moments [Hj], R2, to be less than 
0.90. This ratio is calculated from the momenta of all 
charged and neutral particles in the event. To suppress 
backgrounds from e + e~ — > t + t~ events, which tend to 
have lower multiplicity compared to B° — > 77 events, 
the number of reconstructed charged tracks in the event 
is required to be greater than two. 

We define a signal region in the toes~A£ plane by fit- 
ting each variable in B° —> 77 MC and selecting a range 
around the peak of the distribution. The toes distribu- 
tion is fit with a Crystal Ball (CB) shape |l3[ and we 
take a ±3 a region around the CB peak corresponding to 
5.27 < toes < 5.29 GeV/c 2 . The AE distribution is pa- 
rameterized using a double-sided modified Gaussian with 
tail parameters given by 



f(AE) oc exp 



-{AE-tf 



2o L.R 



a L:R (AE - /i) 2 



(1) 



where /i is the peak of the AE distribution, ctl,r are the 
distribution widths, and ol.r are tail parameters on the 
left and right side of the peak, respectively. We again 
take a ±3 a interval around the peak corresponding to 
-0.30 < AE < 0.13 GeV. This region is blinded in 
the on-resonance data until the maximum likelihood fit 
is performed. With this definition of the signal region, 
an toes sideband region is defined as 5.20 < toes < 5.27 
GeV/c 2 , and lower and upper AE sidebands are defined 
as -0.50 < AE < -0.30 GeV and 0.13 < AE < 0.50 
GeV, respectively. 

We select photon candidates according to criteria cho- 
sen to remove poorly reconstructed photons. The energy 
of the cluster must be spread over at least 10 crystals 
with all crystals in the cluster having active electron- 
ics with correct calibrations. The shape of the cluster 
must be consistent with that of a photon in the defined 
energy range; we therefore require clusters to have lat- 
eral moments 14[ in the range 0.15 < /i a t < 0.50. To 
ensure the shower is fully contained within the EMC 
volume, only photons whose polar angle is in the range 
22.9° < pg < 137.5° are selected. The photons are kept 
for further analysis if they are isolated from all other 
clusters in the event by at least 25 cm. 

If an e + e~ collision that results in a trigger is accom- 
panied by another e + e~ collision nearby in time, EMC 
signals from the out-of-time collision may populate the 
event of interest. Due to the large Bhabha-scattering 
cross section, these "pile-up" events often involve high- 
energy electrons which may produce EMC clusters. Since 
the tracking detectors are sensitive over a narrower time 
window, electron-induced tracks may not point to the 
EMC clusters, causing them to be treated as photon can- 
didates. If their energy is also measured incorrectly, ar- 
tificial B° — > 77 candidates may result. This scenario 
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is effectively rejected by requiring the total event energy 
to be less than 15.0 GeV and the cluster time of each B 
candidate photon to be consistent with the trigger event 
time. 

The dominant source of backgrounds are photons pro- 
duced from high energy ir° and r\ decays in continuum 
events. These events are suppressed using a likelihood 
ratio rejection technique. Each B candidate photon is 
separately combined with all other photons in the event 
and the invariant mass, m 7 y, and energy of the other 
photon, are used to calculate a likelihood ratio given 
by 



U = 



(2) 



In this equation, i is a label for ir° or rj, and P represents 
a two-dimensional probability density function (PDF). 
For each B candidate photon the pairing that gives the 
largest value of the likelihood ratio is assigned. The 
signal PDF, P s ; g , is constructed using simulated B°B° 
events containing a B° — > 77 decay where all B candi- 
date photon pairings are used. The PDF for a 7r° or rj, 
Pi, is constructed from simulated e + e~ — > qq and e + e~ 
— > t + t~ events. The B candidate photon in this case is 
required to be produced from a tt° or rj decay, while the 
other photon daughter is required to be reconstructed 
in the calorimeter. The energy of the other photon and 
the invariant mass of the pair are then used to construct 
the 7T° and rj PDFs. A likelihood ratio near 1.0 (0.0) is 
consistent with the B candidate photon originating from 
a 7T° or -q (signal B). Figure [5] shows the L^a and L v 
distributions for B° — > 77 MC events and for B candi- 
date photons from 7r° and r\ decays in MC continuum 
background events. 

For high energy tt° decays with E„o > 2 GeV, the 
daughter photons may not be separated enough in the 
EMC to be resolved individually. In this case the photon 
clusters are said to be "merged" . A merged ir° can mimic 
a B candidate photon because the cluster will have the 
full energy of the parent ir° and will have no associated 
track. At a given energy, the second moment of the en- 
ergy distribution around the center of the cluster will be 
different for a photon and a merged 7r°. This allows for 
the construction of a quantity called the merged 7r° con- 
sistency based on the energy and second moment of the 
cluster. This variable compares the two inputs against 
known distributions from photons as well as merged 7r° 
decays to estimate the likelihood that the cluster origi- 
nates from either source. The ir° energy range that con- 
tributes B candidate photons in this analysis begins at 
about 2 GeV and extends above 6 GeV. To reduce this 
source of merged 7r° background, we select B candidate 
photons whose merged ir° consistency is compatible with 
that of a photon. 

In the T(45 l ) CM frame, the B mesons are pro- 
duced nearly at rest and subsequently decay isotropi- 
cally, whereas events produced in continuum processes 
are typically collimated in jets along the qq axis. This 
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FIG. 2: (a) 7T° likelihood ratio for B candidate photons 
in simulated signal events (open histogram) and simulated 
continuum background events (shaded histogram) where the 
photon is required to originate from a 7r° decay. (b) 77 
likelihood ratio for B candidate photons in simulated sig- 
nal events (open histogram) and simulated continuum back- 
ground events (shaded histogram) that are required to origi- 
nate from an r\ decay. Events where both B candidate photon 
likelihood ratios are less than 0.84 are selected, as denoted by 
the arrows. 



difference in event shape is exploited to further separate 
signal from continuum background events using a neu- 
ral network (NN) multivariate classifier. The NN utilizes 
19 input variables that characterize event level features 
whose distributions show separation power between sig- 
nal and continuum background events. The inputs in- 
clude the minimum distance between each B candidate 
photon EMC cluster and all other EMC clusters in the 
event, the polar angle of the B candidate momentum 
in the lab frame, the number of neutral particles recon- 
structed, the number of tracks in the event, the miss- 
ing energy and the total transverse momentum. We also 
include quantities that characterize the rest-of-the-event 
(ROE), calculated using all reconstructed particles in the 
event except one or both B candidate photons. These 
are: the polar angle of the event thrust axis when ei- 
ther B candidate photon is removed; the first, second, 



9 



C\J 

o 



CI 

CD 
> 
LU 



10 b 



icr 



10 1 



10 - 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 



0.8 0.9 1 

NN response 



FIG. 3: Output of the neural network on the validation sam- 
ples of simulated signal (open histogram) and background 
(shaded histogram) events. Selected events have a NN re- 
sponse greater than 0.54 as denoted by the arrow. 



and third angular moments of the event when the B can- 
didate photon with the larger energy in the lab frame is 
removed; and the second angular moment with respect to 
the thrust axis of the event when both B candidate pho- 
tons are removed. Additionally, i?2, the total sphericity, 
and the sphericity with both B candidate photons re- 
moved are used; all energies and momenta are calculated 
in the CM frame. We train and validate the NN using 
independent samples of B° — > 77 and continuum MC 
events. The training samples are constructed by first ap- 
plying the photon quality, event pile-up, and merged tt° 
selections to the MC events. The surviving events are 
randomly divided into one set for training and one for 
validation. Each set contains 45,200 events where half 
are B° — > 77 MC events and the other half are contin- 
uum MC events whose composition of e + e~ — » qq and 
e + e~ — > t + t~ events is scaled to the luminosity of the 
on-resonance data for each component. The parameters 
of the NN are tuned to achieve the highest level of back- 
ground rejection while avoiding overtraining the classi- 
fier. The validation sample is then used to verify its 
performance. The NN response is a value between 0.0 
(background-like) and 1.0 (signal- like) . Figure [3] shows 
the NN response for the validation samples for B° —> 77 
and continuum background MC events. 

The selection criteria for L^o, L v , the NN response, 
and the number of tracks are optimized using B° — > 77 
MC events and on-resonance sideband data. The num- 
ber of background events is estimated by extrapolating 
the sideband data into the signal region. The optimiza- 
tion proceeds by iterating over the space of each variable 
individually until we find a maximum of the figure-of- 
merit given by the expression e s j g /(3/2 + y/B)- In 
this equation, e S i g is the efficiency of the event selec- 
tion derived from B° — > 77 MC and B is the number 
of background events. The iterative process continues to 
cycle through all variables until the selection values con- 



verge. We find the optimum values to be: L^o < 0.84, 
L v < 0.84, NN response greater than or equal to 0.54, 
and the number of tracks to be greater than two. The 
optimum selection criteria are found to have an overall 
efficiency of 26.7% on a collection of 1.96 million simu- 
lated B° — > 77 events, while rejecting about 99.9% of 
background events. 



B. Exclusive B Decay Backgrounds 

Backgrounds from B decays that may peak in the Toes 
and AE signal region are studied with large samples of 
simulated events. Twelve B decay modes were identi- 
fied as potential background sources, and exclusive MC 
samples were generated for each mode. The optimized 
event selection is applied to each of these samples and the 
estimated number of background events expected in on- 
resonance data are determined from the latest branching 
fractions [ill- After scaling the yields of these modes to 
the luminosity of the on-resonance data, it is estimated 
that they contribute a total of 1.18 ± 0.22 background 
events to the signal region. This number is comparable to 
the expected number of signal events predicted from the 
SM branching fraction (~ 4 events). The modes expected 
to contribute significantly are B° -> 7r°7r°, B° -> n°r), 
B° — » r)r), and B° — > ory. The toes distributions of 
these modes peak at the same value as true signal events, 
while the AE distributions peak at a value less than zero. 
This difference in shape of the AE distributions between 
B° —> 77 decays and these "peaking" background B de- 
cays is exploited by adding a component that describes 
them to the maximum likelihood function. 



IV. MAXIMUM LIKELIHOOD FIT 

The signal yield is extracted using a two-dimensional 
unbinned extended maximum likelihood (ML) fit in the 
region to es > 5.2 GeV/c 2 and -0.5 < AE < 0.5 GeV. 
The likelihood function for a sample of N events with 
signal, continuum, and peaking BB background compo- 
nents is given by 




(3) 



where i in this equation is an index for the three compo- 
nents in the fit and rii is the event yield for each. Since 
the correlations between Toes and AE are found to be 
small, the signal and continuum background PDFs, Vi, 
are each defined as a product of one-dimensional PDFs in 
the observables Xj € {toes, AE}, with parameters d?j. A 
two-dimensional histogram PDF is used for the peaking 
background component. 

The signal PDF shapes for toes and AE are deter- 
mined from simulated £>° — > 77 events. The toes distri- 
bution is parameterized by a Crystal Ball function [l3j], 



and the AE shape is parameterized by a double-sided 
modified Gaussian with tail parameters given by Equa- 
tion (JlJ. In the ML fit, the signal PDF parameters are 
fixed to the MC-determined values. All fixed signal pa- 
rameters are later varied to evaluate the systematic un- 
certainty that this choice of parameterization has on the 
signal yield. 

The continuum background toes distribution is param- 
eterized by an ARGUS shape [Tg, while the AE distri- 
bution is fit with a first-order polynomial. The endpoint 
of the ARGUS function is fixed to the kinematic limit for 
B decays (5.29 GeV/c 2 ), while all other parameters are 
allowed to float. The PDF for the peaking background 
component is parameterized using large samples of simu- 
lated exclusive B decays in the form of a two-dimensional 
histogram PDF in toes an d AE. Both the shape and 
yield of this component are fixed in the ML fit. The 
yield is fixed to 3.13 ±0.54 events, which is the predicted 
number in the fit region determined from the exclusive 
MC studies. The fixed peaking background PDF shape 
and yield are later varied to evaluate the systematic un- 
certainty on the signal yield. 

The fit is validated on an ensemble of prototype 
datasets whose signal and background content and shape 
are as expected in the on-resonance data. For the signal 
content, one-half, one, and two times the expected SM 
branching fractions [l| are used to populate the prototype 
datasets. Two types of datasets are constructed: one 
where both the signal and background events are gener- 
ated by randomly sampling from their respective PDFs, 
and the other where the background events are generated 
from a random sampling of the background PDF while 
the signal events are embedded directly from the simu- 
lated signal dataset. The results of the validation show 
that the fit returns the correct signal yield and does not 
introduce a bias. 

The on-resonance T(4S) data contains 1679 events af- 
ter the optimized event selection criteria are applied. 
We perform the ML fit to extract the signal yield and 
find A s i g = 2ll^2 events corresponding to a statisti- 
cal significance of 1.9 a. The significance is computed 
as \/2 • Aln£, where Aln£ is the difference in the log- 
likelihood between the best fit to on-resonance data and 
a fit where the signal yield is fixed to zero. Figure 2] shows 
projections of the PDF components from the ML fit. For 
the toes projection, the range of AE has been restricted 
to -0.30 < AE < 0.13 GeV. For the AE projection, the 
range of Toes is restricted to toes > 5.27 GeV/c 2 . 



V. SYSTEMATIC UNCERTAINTIES 

Systematic errors that affect the calculation of the 
branching fraction are investigated and include uncer- 
tainties on the number of B°B° events in the dataset, 
signal efficiency, and the signal yield from the fit. Differ- 
ences between data and MC can lead to an error on the 
derived signal efficiency. The identified sources that can 
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FIG. 4: Projections of the ML fit onto tties and AE. (a) 
The projection of the tubs component when the range of AE 
has been restricted to -0.30 < AE < 0.13 GeV. (b) The 
AE projection of the fit when the range of ?ties has been 
restricted to jites > 5.27 GeV/c 2 . The points represent the 
on-resonance data. The solid curve represents the total PDF, 
the dashed curve is the continuum background component, 
the dot-dashed curve is the signal component, and the long- 
dashed curve is the peaking background component. With 
an expected yield of approximately one event, the peaking 
background component is nearly indistinguishable from the 
x-axis. 



lead to this error include uncertainties in tracking, track 
multiplicity, photon reconstruction, the L v o and L v re- 
quirements, and the NN selection. The uncertainties in 
the modeling of the signal and BB background shapes 
in the maximum likelihood fit can also affect the uncer- 
tainty on the signal yield. 

The systematic error associated with counting the 
number of BB pairs in the dataset is 1.1%. The num- 
ber of B°B° pairs is obtained by multiplying the TjiS) — > 
B°B° branching fraction, equal to (48.4±0.6) % [UJ, with 
the total number of BB pairs, from which a systematic 
error of 1.7% is assigned. A study of the track finding 
inefficiency results in an assignment of a 0.2 % system- 
atic uncertainty for the selection of events with at least 
three reconstructed tracks. The uncertainty in the signal 
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efficiency due to the requirement of at least three recon- 
structed charged tracks is estimated to be 3.4%, includ- 
ing components for both generator-level simulation er- 
rors and detector-associated data versus MC differences. 
The efficiencies in data and MC for detecting high-energy 
photons which pass the selection criteria were compared 
using a sample of e + e~ — > /i + /i~7 events with the pho- 
ton energy in the CM frame restricted to be consistent 
with the energy of B° —> 77 photons. No significant dif- 
ference is observed. An uncertainty of 2% per photon is 
assigned to account for possible data versus MC differ- 
ences due to the required minimum distance between the 
candidate-photon cluster and all other clusters, based on 
a study that embedded high-energy photons in both data 
and MC events. We combine this uncertainty linearly for 
both photons in B° — > 77 and assign an overall photon 
efficiency systematic of 4%. The cluster time selection 
is compared in data and MC and we assign a systematic 
uncertainty of 0.7% for each B candidate. The system- 
atic uncertainty due to the 7r° and 77 likelihood ratios is 
estimated to be 1.0% for each, based on a study that 
embedded signal-like photons in data and MC events in 
which one B meson was fully reconstructed in the de- 
cay B — > Dir. The signal efficiencies for the L^a and 
selections for data and MC are calculated by pairing 
the embedded signal-like photon with all other photons 
in these events that are not associated with the recon- 
structed B. The systematic uncertainty due to the NN 
is estimated by comparing the efficiencies of data to MC 
in signal-like events. Signal-like events are selected by 
applying all event selection criteria, but reversing either 
the NN, the L^o or the L v selection for one of the B 
candidate photons. The efficiencies are then calculated 
from events in the fit region with the signal region ex- 
cluded. For all selection reversal scenarios, the ratio of 
the efficiencies is found to be consistent with unity and 
has a typical statistical error of 3.0%, which is taken as 
the systematic uncertainty. 

The systematic uncertainty on the signal yield due to 
the choice of fit model has four components. The first 
is due to the fixed signal shape for the tttes and AE 
PDFs. To estimate the systematic uncertainty related to 
this choice of parameterizations, the fixed parameters are 
varied within their ±1 a errors, the on-resonance dataset 
is refit, and the change in the signal yield is calculated 
for each parameter. The total systematic uncertainty 
is taken to be the sum in quadrature of all variations 
and is found to be 0.6 events. A comparison between 
photon response in e + e~ — > fi + /i~7 events in data and 
MC, using photons in the energy range relevant to the 
decay B° —> 77, shows that the size of the variation 
in the signal shape parameters is sufficient to take into 
account any systematic effects from parameterizing the 
signal PDF shapes from MC. The second component is 
due to the parameterization choice for the signal shape. 
While the PDF used to fit the 77t.es distribution repli- 
cates the shape in B° —> 77 MC, there is a slight dis- 
agreement between the AE distribution in B° — > 77 MC 



and that of the double-sided modified Gaussian used to 
describe it, Eq. (fT]). To test how large of an effect this 
difference may have on the signal yield, the AE distribu- 
tion is parameterized using a Crystal Ball shape [l3[ that 
provides a larger discrepancy from the MC shape. En- 
sembles of simulated experiments are performed wherein 
each experiment consists of an independent dataset fit 
first using the AE parameterization described in Sec- 
tion IIVI and then using the alternative parameterization 
described here. The signal yields for each fit are com- 
pared and the average difference is found to be 0.2 events 
which is assigned as the systematic uncertainty. The 
third component is due to the choice of the AE contin- 
uum background shape. Repeating the fit to data using 
a second-order polynomial for AE results in an increase 
of 1.9 events in the signal yield, which we take to be the 
systematic error for this component. The fourth compo- 
nent is due to the choice of shape and normalization for 
the peaking background PDF, both of which are fixed in 
the ML fit. The shape is fixed from the tttes and AE 
distributions of simulated exclusive B decays, while the 
yield is fixed to the expected number of events in on- 
resonance data. To estimate the systematic uncertainty 
on the signal yield, both the peaking background yield 
and shape are varied. The uncertainties on the yields 
of the individual peaking background modes are added 
linearly to determine the uncertainty on the total yield 
of the peaking component in the ML fit. This results 
in a range for the peaking yield between 2.02 and 4.24 
events. The shape of the peaking PDF is varied by re- 
placing it solely with the shape derived from B° —> 7r°7r° 
MC. The AE distribution of this mode most closely re- 
sembles that from B° — > 77 MC. Another ensemble of 
simulated experiments were performed where the differ- 
ences in the signal yield were compared when using these 
set of extreme variations in the peaking component and 
those from the ML fit to data. We take the maximal 
change in the signal yield and assign a conservative sys- 
tematic uncertainty of 0.5 events. The four components 
are added in quadrature to give a systematic uncertainty 
for the ML fit of 2.1 events corresponding to 9.9%. 

These results are added in quadrature to give a total 
systematic uncertainty on the signal yield of 2.6 events, 
corresponding to 12.1%. The separate contributions to 
the systematic uncertainty are summarized in Table |TJ 

VI. RESULTS 

The branching fraction is calculated from the measured 
signal yield using 

B(B° -> 77) - ; f S % , (4) 

where N s i g is the signal yield from the maximum likeli- 
hood fit, £ s ; g is the signal selection efficiency determined 
from simulated B° — > 77 events, and N BO g is the num- 
ber of neutral B meson pairs in the on-resonance dataset. 
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TABLE I: Summary of the systematic uncertainties expressed 
as a percent of the signal yield. 



Source Uncertainty on N s i s (%) 



B U B" counting 


1.7 


Tracking efficiency 


0.2 


Track multiplicity 


3.4 


Photon efficiency 


4.0 


Cluster time 


0.7 


L n o and L v 


2.8 


Neural network 


3.0 


Fit uncertainty 


9.9 


Sum in quadrature 


12.1 




FIG. 5: Likelihood curve from the maximum likelihood fit, as 
a function of the signal yield, after convolution with a Gaus- 
sian shape whose width is equal to the total systematic error. 
The shaded region corresponds to the integral of the curve 
up to 90% of its total area. The y-axis, L/Lo, is the ratio 
of the likelihood function for a given JV s i g to the maximum 
likelihood. 

We calculate the branching fraction to be 

6(5° ^77) = (1.7 ±1.1 ±0.2) x 1(T 7 , (5) 

with a statistical significance of 1.9 cr, where the first error 
is statistical and the second is systematic. 

The upper limit at the 90% CL is obtained by inte- 
grating the likelihood curve resulting from the ML fit 
from zero to the value of N a i s which contains 90% of 



the area under the curve. To incorporate the systematic 
uncertainty into the determination of the upper limit, 
the likelihood curve is convolved with a Gaussian shape 
whose width is equal to the total systematic uncertainty 
of 2.6 events. This yields a value of iV S ig = 39 events 
corresponding to an upper limit of 

B(B° -> 77) < 3.3 x 10~ 7 (90% CL). (6) 

Figure [5] shows the likelihood curve from the fit after 
convolution with the Gaussian shape. The shaded region 
corresponds to the 90% integral of the curve. 

This limit is nearly a factor of two below the best pre- 
vious upper limit of B(B° — > 77) < 6.2 x 10 -7 set by 
Belle [J], and is consistent with the SM branching frac- 
tion. This limit may allow tighter constraints on models 
that incorporate physics beyond the SM. 
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